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1. Introduction 

One major application of lattice QCD (LQCD) is the determination of hadronic matrix ele¬ 
ments of electroweak processes. These calculations are necessary for interpreting the results of 
experimental searches for physics beyond the Standard Model as well as an improved understand¬ 
ing of Standard Model predictions. In the last few years, we have seen the emergence of controlled 
calculations of processes involving multiple (two) hadrons in the initial and/or final state. Unlike 
the single hadron ground state spectrum or matrix elements, there is no simple relation between 
the finite-volume (FV) matrix elements and the infinite volume (°oV) transition amplitudes. The 
complications arise from the fact that the multi-particle states can not be isolated asymptotically as 
single-particle non-interacting states, and the finite cubic box does not respect the full 0(3) spatial 
symmetry but only the cubic symmetry, so angular momentum is no a good quantum number. For 
recent reviews of this topic, see Refs. [JT], |2 ]. 

In this talk, we review recent formalism that has been developed to study a wide class of 
1 —>■ 2 processes in the presence of an external current [Q]. We generalize the original work by 
Lellouch and Liischer, which was derived for K —>• Tin decays 0 with specialized kinematics. 
This pioneering work was extended to all elastic channels under the inelastic threshold [j|]. For 
this S-wave kaon decay, in the rest-frame of the kaon, the partial wave mixing induced by the 
periodic cubic box is highly suppressed and was therefore not considered in these first works. We 
were motivated to determine a master formula derived with as few approximations as possible. We 
present this master equation, relevant for pseudo-scalar particles where the final state is composed 
of any number of coupled channels. We also demonstrate how partial-wave mixing is manifested 
in the determination of the infinite volume transition amplitudes from the finite volume matrix 
elements. 1 


2. Two-particles in a box: a review 


We begin with a brief review of the now well known quantization condition for two-particles 
in a finite volume with zero intrinsic spin, any number of open and coupled channels and any total 
momentum. This is a generalization of the well known two-particle Liischer formula [j7[ [sj] which 
was later extended to boosted two-particle systems ©IH- The formula presented here was first 
obtained in Refs. EMI: 2 


det 


K(£„) + (F y (£„))' 


= 0 . 


( 2 . 1 ) 


In this equation, IK is the two-particle K-matrix, which is directly related to the scattering amplitude, 
. via ./£ 1 = IK 1 — /'P 2 /2, where P 2 is the phase space factor, which is a diagonal matrix in the 
open channels and is defined in Ref. [|3|]. Although the K-matrix is diagonal in angular momentum 
space, the finite volume matrix, ¥ v , is not. In general, this is a matrix that depends on the geometric 
and boundary conditions of the volume, as well as the energy and momenta of the system. 3 In a 

1 A similar work determining the volume modifications for yN —> JtN was performed for restricted kinematics and 
no partial wave mixing, but with the inclusion of spin in Ref. [^] 

2 The discussion in this section follows closely that of Ref. S 

'The formula for F v is generalized to systems with arbitrary intrinsic spin, in a volume that is a generic rectangular 
prism with twisted boundary conditions in Ref. |[i~3||. 
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cubic box with periodic boundary conditions 


IX; ] ImrJ'm,, 


8 kP* 




0 ,M l" ,m" K'jiOti 


where the function cf m is defined as 


4.^-f (t 


1—2 


&* m [\-(k*L/2nf 


%L\^A = L 


|r| / L/ m (r) 


rt»S rl - x2 A 


( 2 . 2 ) 


(2.3) 


and kj on is the on-shell relative momenta of the two particles in the center of mass (c.o.m.) frame. 
The determinant in Eq. f p~T| ) is evaluated over the inner product space of angular momentum and 
open two-body channels. 4 

It is useful to define the residue of the fully dressed two-particle propagator in finite volume in 
terms of the matrix M (the argument of the determinant in Eq. ( |2.![ )) 

M(Po,m) = ^(E 0 ,m) + (^(Pq.m)) • (2.4) 


The residue is then given by 

&A f ,n f = K adj [M(Po,m)] tr 


adj [M(P 0 ,m)] 


5M(P 0 „ 


M) 


dP, 


0 ,M 


K. 


(2.5) 


Po,M —j ,/ly 


where adj[M] = M 1 det[M] is the adjugate of the matrix M. E,\ fMf is the ray eigenvalue of the Ay- 
irreducible representation (irrep) of the two-particle system. 


3. 1 —y 2 matrix elements 


With the above formalism in hand, it is straightforward to arrive at the generalized Lellouch- 
Luscher result. Let be the infinite volume transition amplitude for a 1 —» 2 process under 

the presence of an external current, where (A/,«/) denote the quantum numbers discussed above 
for the two-particle state and A c labels the quantum numbers of the external current. 

Define | {E Af ,„ f Pf,L\ ^ (0, P y - P,-) |Pa,oP*;L) | as the finite volume matrix elements of the 
external current. The resulting non-perturbative master equation relating these two quantities is 


I (E Aft n f Pf,L\ (0, P f - Pf) |Pa, 0 P«-;P) I 


2Pa,,o 


A f,n f \Xc 


(3.1) 


where c/tf^/fy) is the normalization of the initial(final) finite volume state, |£a ; ,oP/;P) ( \E\ f ^ nf Yf,L )), 
which in Ref. [|3|] was set to 1. The subscript A ; - denotes the irrep of the single particle groundstate, 
which corresponds to n, = 0. Equation [3. l[ which is the main result of this work, relates the finite 
volume matrix to, in general, an infinite number of transition amplitudes, .s/ is understood as a col¬ 
umn vector (and ,o/ f a row) in the spaced defined above (angular-momentum space © open channel 
space). Just as with the two-particle quantization condition, in practice one is required to performed 
a truncation over the contributing elements of the transition amplitude. For low-energies where the 
higher partial waves are suppressed, this is a justified approximation. For a detailed derivation, see 
Ref. [§. _ 

4 For a very nice application of this coupled channel formalism, see the work of the Hadron Spectrum Collaboration 
computing the scattering parameters in the KK-Krj system [[Id], p5|. 
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3.1 K-^- 7171 

In the original work by Lellouch and Liischer [ 4 ]]. the initial and final states were assumed to 
be at rest. To recover this result, one first truncates the quantization condition, Eq. ( p~T| ) to a system 
of S-wave interactions with suppression to all higher partial waves, 

cot 8 S + cot 0 o„ = 0, (3.2) 


where we have introduced the pseudo-phases tyf defined in terms of the cf functions, Eq. (2.3), 


4% 


<7a cot (j)f m = - —cf m (qf ; L). 

<iL 


(3.3) 


The ratio of the infinite transition amplitude to finite volume matrix element is then approximated 


Ws,. 


n f \ 


\6tzE k E* f d(8 s + 0 


\{7t7t, E nf V , \fjXf\L\ ^weak (0,0) |/T, E K V\L) | 2 q*£ dP 0M 


(3.4) 


Po,M—En 


Ek is the energy of the incoming kaon and the symmetry factor q is equal to 1/2 if the final two 
particles are identical and equal to 1 otherwise. Using Eq. ( |2.2p . one can obtain the limit in which 
the final two particles do not interact, 


«/1 


2 ££_ = 2FilEl l 3 . (3 .5) 


\(7:n,E nf P,A f iif,L\jr weak (0,0)\K,E K P-,L)\2 f ree % 

In the limit of two identical final state pions with E n . = Ek, one recovers Ref. 0. 


3.2 One nontrivial example: tzj* —>• Tin and B Tin + it 

Consider a final state with two degenerate pions in a P-wave. Ignoring contamination from the 
F-wave, the spectrum will satisfy the equation 

cot 8p + cot (f)p = 0, (3.6) 

where the pseudo phase <f)p can be written as a terms of the (j)f m defined in Eq. ( p.3| ) 

cot (f)p = ^COt </>,?„ + a 2 0j A cot 020 + «22,A C Ot 022) • ( 3 - 7 ) 

where a/o a are constants and have been tabulated in Table II of Ref. 0 for irreps that have overlap 
with P-wave for symmetry groups corresponding to boosts v/d 7 < 2. The resulting ratio of infinite 
volume transition amplitude to finite volume matrix element is then 

_ IstXfX? _ _ 1 16 %E X E* f d(8 P + 0p) 

|(a:7r,E„ / P / ,A /j u / ;L| c / Ar (0,P / -P()|7r,£',P;;L)| 2 % q* f dP 0 , M 

where Ay compactly encodes the labels associated with the final state that are explicitly discussed 
in Ref. 0. Unlike the scenario considered before, here the external current can inject arbitrary 
four-momenta and the symmetry factor q has been set to one since the two particles must be 
distinguishable. This result is directly applicable for studying processes such as Tty* -A Tin and 


,(3.8) 
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B —>• nn + M. For a decay process such as the one considered in Ref. [[17[]. B —> nK + 11, the 
final two hadrons are no longer degenerate and therefore S and P-wave will in general mix when 
the system is boosted. For these systems, one can still use Eq. (fO|). See Ref. [|3j] for a detailed 
discussion of irreps and partial wave mixing in this system. 

In this case, the non-interacting limit corresponds to Eq. ( fL5| ) with q = I. For strongly inter¬ 
acting systems with rapidly varying phase shifts, the right hand side of Eq. ( [3.8| ) can largely deviate 
from it the free limit. Near an infinitely narrow resonances, the derivative of the phase shift will 
diverge. Parameterizing the phase shift near 90° as a relativistic Breit-Wigner 


tan 8 


e* r p 

(m 2 — E* 2 ) ’ 



E*—m p 


— sin” 8 cot 8 
dE 


E*=m p r P m P 


(3.9) 


where E p = ^P 2 + m 2 . As the width goes to zero, the derivative divergences, as expected. Ignoring 
contributions from the derivative of the pseudo-phase, which only diverges near the free-particle 
poles, one obtains the following result for the right hand side of Eq. (3.8) near a resonance mass 
with a narrow width 

Wx r x\ 2 


1 1 


\(Tnc,E nf Vf,KfHf,L\/ K {Q,Vf-?i)\n,EiVuL)[< 


2 Er 


1 6nE n 


€ q* P 


(3.10) 


For a system with a resonance, one can parametrized the transition amplitude for aX^ nn, where 
X is stable single particle state and the nn system has been projected onto the p-channel, in the 
following form 


Xf \X C 


= f: 






(£■*, |Q| 


mt 


\JE*Y(E*) 1 

-E 2 m -iE*r(E *)^| 



(3.11) 


Note that although the numerator of the scattering amplitude in the presence of a resonance is pro¬ 
portional to r (E*), the transition amplitude is proportional to > JT(E *). This can be understood by 
the fact that in a transition amplitude, only one of the external legs couples to the two-particle states, 
while in a scattering amplitude both incoming and outgoing states are composed of two-particle 
states. Near the resonance mass, the transition amplitude divergences inversely proportional to the 
square root of the resonance width 


^X P K>IQD 

—i\F>\J q*pEp/%n 


(3.12) 


where F^^(m p , |Q|) is the X —>• p transition amplitude evaluated at E* = m p . 

In this limit, one finds that the finite volume matrix element calculated via lattice QCD is equal 
to the infinite volume X —> p transition amplitude up to the standard normalization of the states, 


(nn,E nf Pf,AffXf-,L\ ^^.(0,P/ Pi)\n,EiPp,L) 


I'vfK.IQI 


2 xEi+Yf 


(3.13) 


2Ej 2E p 

It is important to emphasize that this approximation only holds when one knows that a resonance is 
very narrow (r r/ttir <C 1) and the energy level determined corresponds to the resonances mass, up 
to small correction that scale with the width. To reliably asses the validity of this approximation, 
one must first determine the phase shift as a function of the energy using Eq. (|Q|). 
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3.3 Comment on the calculation of B —> K*£ + £ with staggered fermions 


Promising observables to search for new physics are the rare decays B —>• K*i + £ and B s 


<pi~£ . These processes were recently addressed with lattice QCD calculations Qlg. I_7] utilizing 


staggered fermions. Unfortunately, due to the resonant nature of the final state, the use of staggered 
fermions (or mixed-action/partially-quenched calculations) present a technical challenge which 
likely invalidate the results for direct comparison with physical transition amplitudes. The problem 
is that the final state K* can decay strongly to the two-particle Kn system, including disconnected 
quark-flow diagrams. In a calculational scheme that does not completely respect unitarity, such 
as the rooted staggered formulation (or mixed-action/partially-quenched schemes), these strong 
decays will contain so called hair-pin interactions with unphysical double poles [|T8|,[T9|,|2(i|]. These 
hair-pin interactions in the s-channel diagrams were shown to invalidate the know relation between 
two-particle energy levels in a finite volume and the infinite volume scattering phase shift [pd| |, 
as the optical theorem is not satisfied even below the inelastic threshold. Scattering in maximally 
stretched isospin channels is protected from these issues in the elastic region, but suffers the same 


problem above inelastic thresholds Q22|]. Consequently, as the K* final state in the above mentioned 
calculations | |~Ut| . p~7| ] can strongly decay, we do not know formal relation between the finite-volume 
matrix elements which were computed and the corresponding infinite volume transition amplitudes, 
nor do we know how to quantify the error made in the calculations. For these reasons, processes 
such as these require a numerical formulation which respects unitarity, at least at this time. 


4. Conclusion and discussion 


In this work we have reviewed recent work that allows for the determination of transition 
amplitudes for 1 —>• 2 processes where the individual particles have no intrinsic spin. This is relevant 
for meson photo production and heavy meson decays for example. We have found a universal 
equation relating finite volume matrix elements and infinite volume transition amplitudes, Eq. (fO[). 
We have reviewed two simple examples on its implementation, K —>• nn and X —» nn (where X is 
a stable single particle state and final state is projected onto a P-wave). Furthermore, we discuss 
the free limit and the narrow width limit of our results. 


The steps needed to reliably determine 1 —> 2 transition amplitudes via lattice QCD are: 

1. Determine two-particle spectrum: For systems involving resonances this requires utilizing 
a wide range of operators to maximally disentangle the energy levels. For example, in the nK- 
Kr\ calculation by Wilson et. al [[Mj. [i~5| l, this required using “ qq-like ”, “ nK-like ”, and “ Krj-like ” 
operators which requires the evaluation disconnected diagrams. 

2. Plug the finite volume spectrum into Eq. (|2.7[) to determine the K-matrix 


3. Parametrize K-matrix: By parametrizing the K-matrix one can then analytically take its 
derivative as a function of the energy of the system and consequently determine the two-particle 
residue using Eq. (2.5). 

4. Evaluate finite volume matrix element via lattice QCD 


5. Determine the 1^2 transition amplitude using Eq. f|3. I\. 
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